Magnetic Coulomb Fields of Monopoles in Spin Ice and Their Signatures in the Internal Field 

Distribution 
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Fractionalisation - the breaking up of an apparently indivisible microscopic degree of freedom - is one of 
the most counterintuitive phenomena in many-body physics. Here we study its most fundamental manifestation 
in spin ice, the only known fractionalised magnetic compound in 3D: we directly visualise the l/r^ magnetic 
Coulomb field of monopoles which emerge as the atomic magnetic dipoles fractionalise. We analyse the inter- 
nal magnetic field distribution, relevant for local experimental probes. In particular, we present new zero-field 
NMR measurements which exhibit excellent agreement with the calculated lineshapes, noting that this experi- 
mental technique can in principle measure directly the monopole density in spin ice. The distribution of field 
strengths is captured by a simple analytical form which exhibits a low density of low-field sites — in apparent 
disagreement with reported /iSR results. Counterintuitively, the density of low-field locations decreases as the 
local ferromagnetic correlations imposed by the ice rules weaken. 



Introduction — The magnetic field set up by a spin con- 
figuration is the most direct manifestation of the underlying 
magnetic moments. The discovery of a new spin state thus 
holds the promise of generating - and revealing its existence 
in - novel properties of the field it sets up. 

A case in point is spin ice which - uniquely among 
magnetic materials in three dimensions - exhibits an emergent 
gauge field and magnetic monopole excitations |Z!|, which 
have analogies in magnetic nanoarrays Isf-Sl- The spin ice 
state has the great advantage of exhibiting phenomena of fun- 
damental conceptual importance in a setting which as we de- 
scribe below is simple enough to be easily and intuitively vi- 
sualised: an order of topological nature manifests itself in the 
fractionalisation of the microscopic dipole degrees of free- 
dom, leading to the deconfined magnetic monopoles ff]. 

Neutron scattering experiments, which provide magnetic 
field correlations in reciprocal space, have produced some of 
the strongest evidence so far for the gauge structure IHQ] and 
'Dirac strings' ifioll that emerge at low temperatures. Another 
probe which has prominently been employed is muon spin ro- 
tation ifllUlsIl (/iSR), which like NMR is sensitive to the local 
fields in real space. For such local probes, studies of the level 
of detail characteristic of the neutron analysis are still lacking. 

We remedy this situation by computing the spatially re- 
solved distribution of internal fields in spin ice. Most fun- 
damentally, according to Ref. 2, the internal fields in spin 
ice contain a contribution from the underlying magnetic 
monopoles. Isolating and identifying this contribution is thus 
of great conceptual importance in corroborating the peculiar 
nature of these unique elementary excitations. 

Here we show how to visualise the monopole contribution: 
by measuring the field strength at the considerably-sized mag- 
netic voids of the lattice we find a radially symmetric signal 
(Fig. nil which is well-described by the Coulomb law, oc 
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FIG. 1. Cross sections of the primitive unit cell. Left: Average 
field strength in 2in-2out spin ice configurations. The logarithmic 
colour scale ranges from dark blue (0.137 Tesla, the smallest average 
field strength recorded) to deep red (6 Tesla). Right: Average field 
strength in completely disordered Ising spin configurations. 



with a coefficient that is in good agreement with the theoret- 
ical prediction iQ] ■ Even if measuring the field strength deep 
inside the sample may be beyond the reach of current exper- 
iments, the Coulomb field of a magnetic monopole near the 
sample surface could be accessible to a sufficiently spatially 
resolved measurement. 

To make contact with ^SR and NMR experiments, we com- 
pute the full field distribution in the unit cell (Fig. This 
provides detailed predictions for NMR experiments, with the 
lineshape (Fig. |2]i in excellent agreement with the first zero- 
field NMR measurements, the results of which we report here. 

Our analysis places strong constraints on the /iSR signa- 
tures of the spin ice state. In particular, it seems highly un- 
likely that the signal detected in Ref. is due to muons im- 
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planted in pristine bulk spin ice. (Our results are consistent 
with earlier estimates of the internal field strength at the muon 
site in [HQ.) 
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FIG. 2. Histograms of the magnetic field strength at the 0(1) sites 
obtained from: experimental NMR data nominally at T — 0.1 K 
(blue) and at T = 0.4 K (cyan) - top axis; Monte Carlo simulations 
in equilibrium at T = 0.6 K (green) and equally weighted ensemble 
of 2in-2out spin ice configurations (red) - bottom axis. The exper- 
imental curves have been shifted so as to match the main peak po- 
sition from numerics (see the Methods Section in the Supplemental 
Information); the vertical axis is chosen to set the maxima equal to 
unity. 

In finer detail, we find that the internal field distribution fits 
well to a simple functional form at all temperatures (Fig. [3]). 
Counterintuitively, we find an enhancement of the weak field 
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FIG. 3. Histograms of the field strengths across a uniform cubic 
grid spanning the primitive unit cell in a system of size i = 4 con- 
taining 16L'' — 1024 spins. Red: without monopoles. Magenta 
and blue: with two monopoles (see text). Green: random spin ice 
configuration (i.e., T much larger than any interaction energy scale). 
Black lines: centres of the super tetrahedra and of the rare earth 
tetrahedra (low and high field curves, respectively). Black triangles: 
P{h) ~ + H§f fit to the random spin ice behaviour In- 

set: Spatial distribution of locations of field strength smaller than 
10 mTesla in at least one of the 10,000 statistically independent con- 
figurations sampled. The dimensionless volume fraction of such sites 
is approximately 5 10^^. 

sites as the temperature is lowered. This is surprising as spin 
ice is a ferromagnet - as defined by the sign of its Weiss tem- 
perature - and one might expect enhanced internal fields to 
appear as spins align. We interpret this unusual behaviour as 
a result of the interplay of the nanoscopic lattice structure of 



spin ice and the slow decay of monopolar fields. 

Overall, our analysis plugs two gaps: firstly, the conceptual 
one between the effective long-wavelength emergent gauge 
theory Hf] and the nanoscale physics of the lattice; and sec- 
ondly, the practical one between theory and real-space exper- 
imental probes. 

Distribution of internal field strengths — In spin ice, the 
magnetic dipoles reside on the sites of the pyrochlore lat- 
tice, which consists of corner-sharing tetrahedra. We provide 
details of this structure as Supplemental Information but all 
that is necessary for digesting the following is: (i) in any 
of the exponentially numerous spin ice configurations, two 
spins point into each tetrahedron and two point out; and (ii) an 
(anti-)monopole coiTesponds to a tetrahedron with three spins 
pointing in (out). 

In Fig. [5] we show histograms of the internal field distribu- 
tion P{h) collected across the primitive unit cell for three dif- 
ferent classes of spin configurations (see the Methods Section 
in the Supplemental Information). We consider the cases of 
monopole-free states (red line) and of configurations contain- 
ing two maximally separated monopoles, evaluating the field 
in a primitive cell containing a monopole (blue dots) or half 
way between the pair (magenta crosses). This is compared 
to a random configuration of Ising spins with local [111] easy 
axes, coiTesponding to an infinite temperature state (green). 

In all cases, at small fields, P(h) oc h^, while at large fields, 
P{h) cx h^^. The latter reflects the geometric probability of 
probing the 1 /r^ divergence of h close to a spin. The former 
is a non-trivial result which will play an important role in the 
interpretation of /iSR experiments further below; it implies 
that a site with a vanishing field is not 'special' in the sense 
that even an entirely flat probability distribution for each of 
the three components of the field vector would yield this func- 
tional form for P{h). 

The presence of a monopole is only weakly visible far away 
from it but nearby its effect is felt strongly - statistical weight 
is shifted from low to higher fields, whereas the overall shape 
of the distribution does not vary appreciably. This is highly 
counterintuitive, if one considers that spin ice 2in-2out tetra- 
hedra are ferromagnetically ordered: all the spins point in the 
same direction, to the extent allowed by the local easy axes. 
For a ferromagnet, one would naively expect that the internal 
fields are larger in the 2in-2out arrangement than they are in 
presence of a monopole or otherwise disordered spins. 

There are two reasons why this happens nonetheless. 
Firstly, the characteristic dipolar correlations between 2in- 
2out tetrahedra lead to an unusually large cancellation be- 
tween fields from different tetrahedra. The spins form 'flux 
loops' where the sum of their dipole moments in fact van- 
ishes. These 'flux loops' get broken down as monopoles ap- 
pear, whose field decreases with distance more slowly than 
that of any dipole. This effect is captured by the dumbbell 
model 120, which accounts well for the long-wavelength as- 
pect of the field distribution. [A more detailed explanation 
can be found in a dedicated section in the Supplemental Infor- 
mation for this paper.] 
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However, to reveal the second reason, such a picture needs 
to be supplemented to account for the detailed structure of 
the field distribution on the lattice scale. This exhibits con- 
siderable local structure, as shown in Fig. [T] near the spins, 
and at the centres of tetrahedra, there are large fields in excess 
of 4 Tesla. By contrast, in the voids between the spins, the 
fields average much lower. Most saliently, at the centres of 
supertetrahedra (see Fig.lH inset), the probability of finding a 
low-field site is greatly enhanced (black dots, Fig.[3]l. Indeed, 
the oscillations in this latter curve provide a crucial pointer; at 
these locations, aided by symmetry, the fields of nearby spins 
can cancel locally, leaving a lower characteristic field scale, 
and hence enhanced low-field probability. 

This shows up in the field distribution averaged across a unit 
cell near a monopole (Fig. |3] magenta line), which follows 
the form P{h) ^ h? / [K^ + i^o)^ derived for the distribu- 
tion of fields due to randomly located and oriented spins 11511 . 
Crucially, the value of i^o is reduced compared to that of a 
defect-free configuration (red line). This picture is backed up 
by the good fit of the above equation to a high-temperature 
state corresponding to a collection of randomly oriented [111]- 
easy-axis dipoles (green line), and hence a high density of ran- 
domly distributed monopoles. 

Finally, Fig. [T] directly demonstrates that, along with the 
breaking of ice rules, the spontaneous spatial organisation of 
fields strengths into high- and low-field locations within the 
unit cell is suppressed. 

Average field due to magnetic monopoles — Having anal- 
ysed the spatial distribution of fields inside the unit cell, we 
next turn to visualising the field set up by a monopole. Re- 
call that magnetic monopoles experience a relative force of 
Coulombic nature. We ask: can one also measure the corre- 
sponding magnetic field (/io/47r)((7/r^)? This is difficult for 
two reasons. Firstly, the internal field away from the lattice 
sites varies tremendously between configurations. Secondly, 
the 'Dirac string' emanating from the monopole carries a 
magnetisation, M, which cancels off the field, H, to give a 
net V • -B = /io(V • -ff + V • M) = 0. These two issues 
can be taken care of by (i) averaging over many configura- 
tions (or, in experiment, over time whilst keeping the observed 
monopole position fixed) and (ii) measuring the field at points 
as far away as possible from any lattice sites. (This also min- 
imises the strong near-field of the spins.) 

In Fig. |4] we display the direction of the average fields at 
the centres of the supertetrahedra set up by two stationary 
monopoles, which visually reproduce the expected hedgehog- 
like monopolar field pattern. We have verified that by sub- 
tracting the analytical Ewald-summed field for two point 
monopoles with charge given in Ref. l^, the residual fields 
appear randomly oriented with average strength tenfold sup- 
pressed. 

To be more quantitative, in Fig. |4] we show the average 
field (over 10^ configurations) evaluated along a line joining 
the two monopoles halfway across a system of 128, 000 spins 
(L = 20), along the [001] direction. The Coulomb field pre- 
dictions are borne out, but masked in part by the line passing 




distance from one monopole (units of ^nn^ 

FIG. 4. Top: Illustration of the magnetic field due to the monopoles 
(red and blue spheres) at the centres of the super-tetrahedra of the 
pyrochlore lattice, visualised by unit vectors in the local field di- 
rection (red-blue arrows), for 1024 spins with periodic b.c, aver- 
aged over 10,000 independent configurations. As shown, each super- 
tetrahedron is formed by 4 regular tetrahedra in the lattice. The 
super-tetrahedra centres are (locally) the farthest points from any spin 
on a pyrochlore lattice site. Bottom: Averaged fields along the [001] 
line joining two monopoles (connected blue dots). The leading be- 
haviour is captured, to within 20% error, by the Ewald-summed field 
from two point magnetic charges at the locations of the monopoles, 
with charge from the theoretical prediction 2/i/ad ~ 4.6 /is/A 
(black line). The periodic deviations from the Coulomb form are 
due to spins which sit very close to the line - this contribution is ex- 
plicitly shown in magenta for the spin at the midpoint between the 
monopoles). 



very close to spins, which contribute a strong and periodic de- 
viation (magenta line) from the theoretical curve at positions 
2\/2n, n = 0, . . . , 10 (in units of run). 

It is tempting to speculate that, if a quantum spin ice ma- 
terial were to be discovered where monopole motion can be 
made slower than quantum dynamics of spin-rearrangements 
preserving the ice rules, the internal fields would be due al- 
most entirely to the emergent magnetic monopoles! 

Experiment I: Zero-field NMR at the 0(1) sites — Given the 
ferromagnetic interactions, the centres of the rare earth tetra- 
hedra experience a large internal field of several Tesla (from 
numerical simulations using Ewald-summed interactions and 
fields - see the Methods Section in the Supplemental Informa- 
tion for further details), which as we show next can be used 
to distinguish between tetrahedra which host monopoles, and 
those which do not. 

The centres of the rare earth tetrahedra in the pyrochlore 
lattice are occupied by oxygen ions (customarily referred to 
as 0(1) oxygens). The ^^O isotopes are NMR active and can 
be used to detect the internal fields in a zero field NMR mea- 
surement. Indeed, two of the present co-authors have carried 
out the first NMR experiments on the 'monopole-free' line as 
shown in Fig.|2] Also shown is a comparison to Monte Carlo 
simulations in thermal equilibrium. The agreement of the line- 
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shape is remarkably good. The experimental spectrum is ac- 
tually peaked at 3.4 Tesla, i.e., about 25% below the calcu- 
lated peak. This shift is most likely due to spatial distribution 
of Dy-4f electrons causing deviations from a point dipole ap- 
proximation (multipolar effects) and/or effects of Dy-O chem- 
ical bonding. 

A promising aspect of these NMR measurements lies in the 
possibility of direct detection of monopoles. Indeed, when a 
tetrahedron hosts a monopole, the field at its 0(1) site drops 
by approximately 13%. This effect is 3 times larger than the 
linewidth due to variations in the field at the 0(1) site because 
of neighbouring monopoles or more distant spins. There- 
fore, the relative intensity of the NMR signal at this pair of 
field values provides a quantitative measure of the density of 
monopoles. In practice, small densities may be hard to detect 
above the background. Also, monopoles must not move over 
the time scale scale of NMR spin-echo experiments (a few 
tens of microseconds) to be detected as a distinct resonance 
line. Given the insights from modelling AC susceptibility re- 
sults, which suggest a hopping rate in the range of millisec- 
onds iflH 171. this condition seems comfortably achievable. 



Experiment U: ^SR — The other major probe of the local 
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magnetic fields is /iSR. Indeed, Ref. [IJ has reported using 
/iSR experiments to measure the monopole charge via an in- 
genious analogy to the Wien effect familiar from electrolytes. 
This is currently a contested experiment 11311 and we would 
like to make a few observations on it from the perspective of 
the current paper. 

First, we note that one important feature of the data pre- 
sented in Ref. [ij is that their signal was extracted from sites 
where the muons experience very low fields, of order a few 
milliTesla. This follows from the observation that doubling a 
transverse field of strength 1 mTesla doubles the /iSR preces- 
sion frequency. From our work it is clear that there is a very 
low density of such sites (see Fig. |3]l, located preferentially 
near the centres of supertetrahedra (Fig.[Tli. It thus seems very 
unlikely that the substantial muon signal emanates from sites 
in the pristine bulk. 



Second, the analysis of Ref. [12| invokes two distinct phe- 
nomena: first, the applied magnetic field induces an increase 
in the monopole density; and second, this increase leads to 
an enhanced depolarization rate. The first feature is indeed 
an expected phenomenon for the steady state of electrolytes 
in a field, but for spin ice, it can only occur as a transient: 
since monopole motion magnetises the sample, there can be 
no steady state with a nonzero monopole current. Moreover, 
in thermal equilibrium, the density of monopoles is believed 
to decrease in an applied field. About this our present work 
has nothing to say. 

However regarding the second step, our results indicate 
that a field-induced regime with heightened monopole den- 
sity would be accompanied by a depletion of low-field sites. 
At least this aspect is qualitatively in keeping with the finding 
in Ref. [l2l that the rate of decay of the /iSR asymmetry gets 
larger in the Wien setting as the applied field is increased. 

Third, combining the above observations leads us to con- 



sider the interesting possibility that the signal arises from 
the action of an enhanced monopole density on muons im- 
planted outside the sample. The idea here is that outside the 
sample the monopole fields would dominate over the much 
smaller fields present in their absence — given that the field of 
a monopole decays less slowly than that of an isolated spin. 
[Interestingly, it was very recently suggested in Ref. [i3| that 
the signal from muons inside the sample is lost altogether be- 
cause of large (compared to 1 mTesla) and fast magnetic fluc- 
tuations and that the only measurable signal due to spin ice 
comes from muons implanted outside, sensitive to stray fields, 
which are analysed in Ref. |T^1 

A very rough estimate based on a monopole liquid sub- 
ject to Debye screening lfl9ll suggests that in the tempera- 
ture range T ~ 0.2 — 0.5 K, the magnetic field set up by a 
monopole measured a distance roughly 100 ± 50 A from the 
sample surface both lies in the range relevant for pSR (be- 
tween 0.1 and 1 mTesla) and dominates that set up by an in- 
dividual spin in the sample. Note that - unlike stray fields set 
up by the magnetisation induced by a uniform external field 
- the monopole density grows with temperature, thus provid- 
ing a qualitative discriminant between the two. A combined 
study of temperature-dependent (uniform) susceptibility and 
/iSR experiments therefore looks like the most prominent di- 
rection to make progress on this issue. 

All such considerations point to the need for more detailed 
studies of what happens near the surface of a sample ll20ll - 
e.g., what happens to the surface monopole density (as a func- 
tion of time and field)? Moreover, little is known at present 
concerning the surface of spin ice samples, e.g., what the na- 
ture of the local crystal fields is or how the magnetic lattice 
terminates. It is also worth bearing in mind that equilibration 
in spin ice may be incomplete at low temperatures, as the equi- 
librium monopole density vanishes exponentially - at 70 mK, 
the lowest temperature accessed so far in /iSR measurements, 
their number density is estimated to be such that there is only 
one pair for a macroscopic volume of a quarter of a cubic me- 
tre! 

Finally, we would be remiss if we did not note that our con- 
siderations here have entirely to do with spatially fluctuating 
static fields and their dephasing of muon precession, whereas 
Ref. 1 1 3l proposes that large dynamical fluctuations are present. 
The tension between their proposal and the magnetic suscep- 
tibility data that are consistent with much slower dynamics for 
the spins will need to be resolved before a consistent account 
of spin dynamics in spin ice can be given. 

Conclusions — We have provided a way of visualising the 
magnetic field of a magnetic monopole inside spin ice, yield- 
ing the first nanoscopic real space picture of this fractionalised 
excitation. Our work attests to the reality of the monopolar 
magnetic field not only at long wavelengths but also on the 
lattice scale. 

As experimental proofs of the existence of monopoles move 
in the direction from thermodynamics towards increasingly 
microscopic 'single-monopole' detection, we hope that this 
work will lay the theoretical groundwork for future searches. 
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such as the ones using NMR or local field probes outside the 
sample which we have outlined above. 
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hedra (namely, the local [111] axis, see Fig. |5]l. The combi- 
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SUPPLEMENTAL INFORMATION 
Spin ice basics 

Here, we briefly collate some basic facts about spin ice - the 
reader interested in a more detailed and complete overview is 
referred to Ref . [l| and the forthcoming Ref . 

The canonical spin ice materials are Dy2Ti207 and 
Ho2Ti207. The magnetic ions Dy'^+ and Ho"^"*" live on a lat- 
tice of corner sharing tetrahedra (pyrochlore lattice) and they 
have large effective spins: J = 15/2 and J = 8, respec- 
tively. A strong local easy axis anisotropy forces the spins 
to lie along the lines connecting the centres of adjacent tetra- 



FIG. 5. Left panel: Illustration of the pyrochlore lattice. Right 
panel: Conventional cubic unit cell in spin ice, containing 16 spins 
(dark gray spheres). The side of the cell equals 2\/2''nn, where r„n is 
the pyrochlore lattice nearests-neighbour distance. A primitive unit 
cell (highlighted by thick blue lines) is 4 times smaller than the cubic 
unit cell and contains only 4 spins. The black lines are the bonds of 
the diamond lattice formed by the centres of the tetrahedra and they 
identify the local [111] axes of the rare earth magnetic moments. The 
distance between the centres of two neig hbou ring tetrahedra is given 
by the diamond lattice constant, = ^3/2 r^n- 

nation of large electronic spins and strong easy axes leads to 
large magnetic moments of the order of 10 Bohr magnetons 
per spin. The long range dipolar interactions are dominant. 
To a good approximation, the full Hamiltonian of these sys- 
tems can be written as 



47r 



E 



Si ■ Sj 3(Si • rjj)(Sj • Vij) 



(1) 



where the summation is over the indices i,j of the spins on the 
lattice, Jcx > is the exchange coupling constant, /io is the 
vacuum permeability, and /i ~ lO/is is the magnetic dipole 
moment of a spin. 

The interplay between the strong local easy axes and the 
dipolar interactions favours nearest-neighbour 'pseudospin 
antiferromagnetism' : a spin pointing into a given tetrahedron 
wishes its neighbors on the same tetrahedron to point out and 
vice versa. The resulting ground-state configuration are de- 
scribed by the ice rules, which stipulate that two spins point 
in and two point out of each tetrahedron. Exponentially many 
configurations satisfy the ice rules, leading to an extensive 
low-temperature entropy of S'p « ^ log |, as verified experi- 
mentally by Ramirez et al. [Qj. Hence, despite a ferromagnetic 
Curie-Weiss temperature « 2K, the spin ice compounds fail 
to develop long-range spin order down to a temperature four 
times lower. 

This peculiar degeneracy has two remarkable conse- 
quences. Most basically, the spin correlators are characterised 
by an emergent gauge symmetry, which results from an ef- 
fective conservation law for the magnetisation density which 
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is implied by the ice rules-the effective theory of spin ice is 
an emergent magnetostatics. In such a theory, the elementary 
excitations in spin ice do not take the conventional form of 
extended domain walls. Rather, flipping a spin in a 2in-2out 
state generates two defective tetrahedra which have 3in-lout 
and 3out-lin spins, respectively. These defects are then free to 
move through the system by means of spin flip processes that 
do not generate further defects: the initial dipolar perturba- 
tion (a spin flip) fractionalises into two point like excitations 
(the two defective tetrahedra). Indeed, one can formally show 
that these excitations are truly deconfined fsfl, being subjected 
only to a Coulombic (1/r) interaction potential. 

These defects can be shown to carry a net magnetic charge, 
so that spin ice is the first experimentally accessible setting 
where free point-like magnetic charges can be observed and 
manipulated. 

Flux loops and temperature dependence of the internal fields 

In order to understand this counterintuitive temperature de- 
pendence of the distribution of low field sites in spin ice mate- 
rials, it is useful to note that this system is very unusual in that 
it is a magnet in which/erromagnetism appears in a frustrated 
way. 

(1) Ferromagnetism: at a single tetrahedron level, the dipo- 
lar interactions are minimised by a ferromagnetic arrangement 
of the spins. The well-known 2in-2out ice rules maximise in- 
deed the magnetic moment of a tetrahedron. (2) Frustration: 
the degeneracy resulting from the frustration-induced under- 
determination leads to longer range correlations, where the 
local moments do not align, like in a usual ferromagnet, and 
hence do not produce a macroscopic moment. Rather, the (ef- 
fective) emergent gauge structure is one which suppresses the 
coarse-grained moments locally, [jstl Pictorially and intuitively 
(if not rigorously), this can be understood by the dipoles align- 
ing head-to-tail to form strings which have a preference for 
closing back on themselves. Such strings exhibit zero net 
mono- and dipole moments in a continuum treatment, which 
in turn leads to an effective suppression of the fields they pro- 
duce at large distances. 

When monopoles are present, the formation of the closed 
loops breaks down as monopoles introduce endpoints of such 
loops. [An open loop terminated by monopoles is commonly 
referred to as a 'Dirac string' in the literature.] Therefore, 
the rapid decay of the loop field due to the vanishing of the 
monopole and dipole moments is replaced by the slow decay 
of the Coulomb field. 

The fields due to a distribution of monopoles thus combine 
to give an enhanced typical field strength, and hence a con- 
comitant reduction in the density of low-field sites. This ex- 
plains the observed effect. 

Indeed, this behaviour is immediately apparent in the 
dumbbell model introduced in Ref. 3|. In that model, each spin 
is replaced by a dumbbell of charges located at the opposite 
ends of the lattice vector separating the centres of two adjacent 



tetrahedra. When the ice rules are satisfied throughout the sys- 
tem, two positive and two negative charges overlap at the cen- 
tre of each tetrahedron. This leads to a charge neutral system 
that gives rise to vanishing internal fields. Adding monopoles 
replaces the effectively charge-neutral field-free system with 
one in which charges are distributed throughout, leading to a 
complex superposition of 'stray fields' of all the monopoles at 
any given site of the lattice. 

These arguments neglect features of the short-distance be- 
haviour, which as we show in this work are themselves 
non-trivial. In particular, lattice-scale variations of the field 
are substantial, ensuring the (energetically required) decrease 
of the effective field strength at the sites of the spins as 
monopoles are created with increasing temperature. 

Now, low-field sites occur due to the cancellation of the 
strong fields of nearby spins. The probability of such local 
cancellation is not much influenced by the presence of a dis- 
tant monopole; however, even if such local cancellation takes 
place, it is then spoilt by the above mentioned monopolar 
'stray fields'. Hence, the occurrence of low-field sites is sup- 
pressed by the creation of monopoles. 

Methods 

Statistically independent spin configurations were obtained 
by means of loop updates that do not introduce violations to 
the 2in-2out ice rules and do not alter the location or charge 
of existing monopoles in the system. The internal fields for 
each configuration were obtained using the Ewald summation 
technique lil], with the Dy2Ti207 parameters given in Ref.0 

The average field strengths in Fig. 1 in the manuscript and 
the histograms in Fig. 3 in the manuscript were obtained from 
a uniform cubic grid of 16, 000 points (corresponding to a 
point spacing of < 0.25 A) across a primitive unit cell of the 
system (Fig.|5]l, averaged over ^ 10"* realisations. 

The Monte Carlo simulations that produced the distribution 
in Fig. 2 in the manuscript used single spin flip updates and 
Ewald-summed dipolar interactions between the spins, for a 
system of 3, 456 spins {L = 6). 

The O-NMR experiments requires ^^O, which is the only 
oxygen isotope with a non-zero nuclear magnetic moment. 
We have enriched ^^O content in a crystal of Dy2Ti207 grown 
by a floating-zone image furnace from the natural abundance 
of 0.01% to approximately 50% by annealing the crystal un- 
der 65% "O2 gas at 1000°C for several days. Only the 0(1) 
sites at the 8a positions located at the centers of the Dy- 
tetrahedra experience large internal fields > 3 Tesla, while 
the internal fields at the 0(2) sites at the 48f positions are 
much smaller and do not contribute to the observed signal. 
Moreover the electric field gradient is zero at the 0(1) sites. 
Therefore, the NMR spectrum at the 0(1) site represents the 
distribution of the local fields (directly obtained by dividing 
the frequency by 5.77185 Tesla/MHz, which is the gyromag- 
netic ratio of ^^O nuclei). To construct a broad-band spec- 
trum, each Fourier spectrum of a pulsed NMR echo is summed 
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while the frequency is being swept (Fourier-step-summing 
technique). Low-temperature experiments have been carried 
out by putting the sample into a one-shot liquid ■^He fridge or 
into a dilution fridge using a He mixture. 



Further experiments 

Surface experiments 



field projected onto the spin direction on a lattice site is known 
to be near constant for an 'ideal' interaction [S. V. Isakov, 
R. Moessner, and S. L. Sondhi [§]; for the conventional dipo- 
lar Hamiltonian, we find a broadening to 0.80 ± 0.03 Tesla. 
This field sets the energy scale for creating bound monopole 
pairs at low temperature. This quantity - as well as the 
Coulomb interaction between two monopoles as they sepa- 
rate - is important for triggering the avalanches which have 
received much attention of late 1 13- 12 ] . 



Ideally, we would like to measure the monopole field (Fig. 4 
in the manuscript) by placing probes at the centres of su- 
pertetrahedra. Interestingly, these coincide with the centres of 
tetrahedra of titanium ions. Placing an NMR-active ion there 
would in principle allow access to the Coulombic form of the 
internal fields. 

However, even if this remains a thought experiment for the 
time being, it is clear that a promising strategy is to suppress 
short-distance fluctuations by moving away from spin loca- 
tions even on the scale of a lattice constant (Fig. 4 in the 
manuscript). The obvious way to do this is to probe the fields 
close to the surface of a sample, as discussed in the previous 
section. The spins and strings being confined to the sample, 
one can thus hope to see the farther-ranging monopole mag- 
netic field undisturbed. Indeed, it is such a set-up that was 
proposed by Zhang et al. ^ to measure the image magnetic 
monopole induced by an electric charge near the surface of a 
strong topological insulator. 



Magnetic avalanches 

Our work also sheds light on a range of non-equilibrium 
phenomena recently measured in spin ice 



prominent analogues in artificial spin ice 111 



|1C|, and their 
ll. While the 



field strength throughout the unit cell fluctuates wildly, the 
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